The Canada-UK Deep Submillimeter Survey V: The 
SubmiUimeter Properties of Lyman-Break Galaxies 



(N 
O 
O 
(N 

CD 

in 

(N 
> 

00 



O 
(N 
O 

:^ 

o 
U 



X 



T.M. Wehh\ S. Eales^, S. Foucaud^, S.J. Lilly^ H. McCracken^, K. Adelberger^, C. 
Steidel^, A. Shapley^, D.L. Clements^, L. Dunne^, O. Le Fevre^, M. Brodwin^ W. Gear^ 



ABSTRACT 

We have used SSO^m maps obtained as part of the Canada-UK Deep SubmiUimeter Survey 
(CUDSS) to investigate the sub-mm properties of Lyman-break galaxies (LBGs). We used three 
samples of Lyman-break galaxies: two from the Canada-France Deep Fields (CFDF) survey 
covering CUDSS-14 and CUDSS-3, and one from Steidel and collaborators also covering CUDSS- 
14. We measure a mean flux from both CFDF LBG samples at a level of ^2a of 0.414 ± 0.263 
mJy for CUDSS-03 and 0.382 ± 0.206 mJy for CUDSS-14, but the Steidel et al. sample is 
consistent with zero flux. From this we place upper limits on the Lyman-break contribution to 
the 850^771 background of '^20%. We have also measured the cross-clustering between the LBGs 
and SCUBA sources. From this measurement we infer a large clustering amplitude of = 11.5 
± 3.0 ± 3.0 h~^Mpc for the Steidel et al sample (where the first error is statistical and the second 
systematic), 4.5 ± 7.0 ± 5.0 h-^Mpc for CFDF-14 and = 7.5 ± 7.0 ± 5.0 h'^Mpc for 
CFDF-3. The Steidel et al sample, for which we have most only significant detection of clustering 
is also the largest of the three samples and has spectroscopically confirmed redshifts 

Subject headings: galaxies:formation — galaxies:evolution — galaxies:starburst — 
cosmology lobservations — (ISM:)dust, extinction — submillimeter 



1. Introduction 

Recent work at optical and submillimeter (sub- 
mm) wavelengths has granted us unprecedented 
access to the high-redshift universe. In particu- 
lar, the Lyman-break selection technique (Steidel 
& Hamilton 1993) has provided thousands of star 
forming galaxies at redshifts z^3 and z~4. These 
data, and the results of other optical surveys (Lilly 
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et al. 1996; Madau et al. 1998; Hogg et al. 1998), 
have shown that the global star formation rate (in- 
ferred from the UV luminosity density at different 
redshifts), increases with redshift to Z'^1 and may 
remain constant to at least z~4 (Sawicki, Lin, & 
Yee 1997; Steidel et al. 1999), implying the be- 
ginning of the epoch of galaxy formation occurred 
at z >> 1. 

Studies of the spatial distribution of Lyman- 
break galaxies show that they are highly clus- 
tered even at these early redshifts (Giavalisco et 
al. 1998). This was initially unexpected since the 
clustering of galaxies had been shown to decrease 
with redshift to z^^l (Carlberg et al. 1997; Le 
Fevre et al. 1996) in line with theoretical predic- 
tions, but the strong clustering of LBGs is actually 
a natural consequence of the effects of bias. Kaiser 
(1984) showed that the high peaks of the density 
distribution in the early universe will have been 
highly clustered, and so objects that form from 
these high peaks, clusters at a redshift of zero or 
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galaxies at a redshift of 3, should also be highly 
clustered. 

The high star formation rates (50-100 M© yr^^) 
and comoving density of the Lyman-break galaxies 
make them attractive progenitors of present-day 
elliptical galaxies (Pettini et al. 1998) though their 
masses are still highly uncertain (Sawicki & Yee 
1998; Somerville, Primack, & Faber 2001). How- 
ever, the Submillimeter Common-User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Tele- 
scope ( JCMT) has revealed a population of dusty 
galaxies with implied star formation rates of > 300 
Mq yr^^ (Small, Ivison & Blain 1997; Hughes 
et al. 1998; Bales et al. 1999). The redshifts of 
these objects are still highly uncertain, with es- 
timates of the median redshift of the population 
lying between 2 and 3 (Bales et al. 2000; Barger 
et al. 1999; Small et al. 2000; Yun & Carilli 
2002; Fox et al. 2002). They have similar spec- 
tral energy distributions to today's ultraluminous 
infrared galaxies (ULIRGs) and often show dis- 
turbed morphology or multiple components, im- 
plying they may be the result of galaxy mergers 
(Lilly et al. 1999; Ivison et al. 2000). Both LBGs 
and SCUBA sources have sufficient star formation 
rates to form present-day elliptical galaxies but the 
extremely high star formation rates of the latter 
mean this can be done on the order of 10^ years, 
as the homogeneous properties of local ellipticals 
indicate is the case. 

The nature of the relationship between these 
two populations remains unclear. An obvious sce- 
nario is one in which they form a continuum of ob- 
jects with the bright sub-mm selected sources rep- 
resenting the highest star forming Lyman-break 
galaxies. Adelberger & Steidel (2000) claim that 
by assuming an Lboi^dust/ Ljjv typical of normal 
starbursts, the LBG population can produce the 
bulk of the 850 /xm background. In this picture 
the two populations are the same objects and a 
separate population of highly obscured ULIRG- 
like objects is not needed. The ratio of optical to 
sub-mm emission for the brighter SCUBA sources 
is, however, much less than for the starbursts con- 
sidered by Adelberger & Steidel and so these ob- 
jects almost certainly represent a separate popu- 
lation (Gear et al. 2000; Bales et al. 2000; Downes 
et al. 1999), but it is unclear whether the fainter 
SCUBA sources, with Sgso < 3 mJy, overlap with 
the LBG population. 



Various optical techniques have been used to 
infer the dust content of LBG's. Pettini et al. 

(2001) have used optical line ratios and Shap- 
ley et al. (2001) fitted the predictions of star- 
formation models to optical and near-IR photom- 
etry. Both have concluded that the most intrinsi- 
cally luminous LBGs, which have higher star for- 
mation rates, contain more dust. However, a more 
reliable way of measuring the dust content is di- 
rectly through sub-mm photometry. Chapman et 
al. (2000) have used SCUBA to observe a sample of 
high-SFR LBGs and they estimate that the 850/im 
flux density is at least two times lower than pre- 
dicted from UV colours. Using the sub-mm map 
of the Hubble Deep Field (HDF), Peacock et al. 
(2000) statistically detected the sub-mm flux of 
galaxies with high UV luminosities, and thus high 
star formation rates. They detect a higher mean 
flux of 0.2 ±0.04 mJy (for galaxies with an inferred 
star formation rate (SFR) of 1 /f^M© yr~^). 

This paper examines the relationship between 
Lyman-break galaxies and SCUBA sources in two 
Canada-UK Deep SubmilUmeter Survey (CUDSS) 
fields and is organised as follows. §2 describes the 
sub-mm and optical/C/U data. In §3 we investi- 
gate the sub-mm fluxes of LBGs. §4 discusses the 
dust properties of LBGs that can be inferred from 
the results in §3. In §5 the correlation function 
between the two populations is presented. In §6 
we discuss the results and their implications. 

2. The Data 

We have mapped two areas within the Canada- 
France Redshift Survey fields (CFRS), CFRS03-h00 
and CFRS14-h52, using SCUBA at 850Atm. This 
study was designed to be a blank-field survey for 
sub-mm-selected sources above 3 mJy but also 
produces statistical information for objects below 
this flux. Each map is roughly 6x8 arcmin^ and 
is a mosaic of single jiggle maps. The sub-mm 
data and the goals and results of the CUDSS are 
discussed in detail in Eales et al. (1999, 2000); 
Lilly et al. (1999); Gear et al. (2000); Webb et al. 

(2002) ; Webb et al, in preparation; Clements et 
al., in preparation. 

Surveys for Lyman-break galaxies have been 
performed over both these areas. In the 14-hour 
field (CFRS14+52) we have access to Lyman- 
break data from two separate sources, the Canada- 
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France Deep Field survey (CFDF) (McCracken 
at al. (2001); Foucaud et al., in preparation) 
and Charles Steidel and collaborators (Steidel, C. 
2001, private communication). The Steidel et al 
list contains 86 galaxies within the SCUBA map 
area and most of these have spectroscopically con- 
firmed redshifts. It is an 7?,- limited survey with 
n < 25.5 (see Steidel & Hamilton (1993) for dis- 
cussion of the filter system) and is roughly 0.5 
magnitudes deeper than the CFDF data. 

The CFDF survey, and in particular the con- 
struction of catalogues and limiting magnitudes, 
is described in detail in McCracken et al. (2001). 
Details of the CFDF Lyman-break selection tech- 
nique are found in Foucaud et al. (in preparation) . 
The survey is /-limited with Iab < 24.5 and has 
selected galaxies over the redshift range 2.9 < z < 
3.5. An additional constraint of {V — I)ab < 
1.0 has been introduced which reduces contami- 
nation by stars and elliptical galaxies at 2;~1.5. 
The CFDF team estimate (from simulations) the 
contamination due to stars to be at the 5% level 
and that due to galaxies below z~2.9 to be 15%. 
The total level of contamination of 20% is compa- 
rable to that of the Steidel sample. In the redshift 
range of 2.9 < 2; < 3.5 and with Iab < 24.5 the 
CFDF method recovers 70% of the Steidel et al 
catalogue. Of the entire CFDF survey a subset of 
26 galaxies fall within our 14-hour SCUBA area 
and 29 within the 3-hour field area. 

3. Submillimeter Flux of Lyman Brecik 
Galaxies 

3.1. Statistical Results 

The CUDSS's large, contiguous maps provide a 
unique opportunity to statistically study the sub- 
mm flux of a relatively large number of Lyman- 
break galaxies. Our method is simple: we mea- 
sure the sub-mm flux at each LBG position in the 
SCUBA maps and take a weighted mean of these 
data, obtaining the weights from the noise maps 
for each field (see Eales et al. (2000); Webb et al. 
(2002)). For high signal-to-noise measurements 
one would measure the flux at the peak of the 
beam (or the point-spread-function) but as these 
objects are all below the noise level we cannot lo- 
cate the peak and must use the value at the LBG 
location (which will be, on average, offset from the 
peak as discussed below). As we will show, it is 



quite crucial to first remove all sources from the 
maps that are definitely not associated with an 
LBG. 

Simulations in an earlier paper (Eales et al. 
2000) of the sub-mm data have shown that the 
offset in our maps between the actual position 
of an object and its recovered position (for the 
Ss50nm > 3 mJy) is usually within 6 arcseconds 
with the peak of the distribution occurring at ap- 
proximately 3 arcseconds (also see Hogg (2001) 
for discussion of offsets). The recovered position 
lies further than 8 arcseconds in only 5-10% of 
cases, and so we removed all SCUBA somces from 
the maps for which there was no LBG within 8 
arcseconds. 

For the 3-hour field we removed all of the 

SCUBA sources, for the Steidel 14-hour sample, 
all but three sources and for the CFDF- 14, all 
but one source. These four sources are listed in 
Table 1 along with their identification probability. 
The source removal, or cleaning, discussed in Eales 
et al. (2000), includes the removal of the entire 
beam template. That is, not only is the positive 
fiux from a source removed but so is the negative 
flux due to the chop. This step is vital since the 
beam profile, convolved with a point source, essen- 
tially extents over > 75 arcseconds in RA and can 
therefore interfere in flux measurements of nearby 
sources if not removed. 

We determined the probabilities that these four 
associations are simply chance coincidences in the 
following way. If an LBG lies d arcscc from a 
SCUBA source, the probability that it is unrelated 
to the SCUBA source is P = 1 — exp(—wn(f), 
where n is the surface density of LBGs (the more 
sophisticated analysis which takes account of the 
the magnitudes of galaxies is described in Webb 
et al. (in preparation). The probabilities, which 
are given in Table 2, are not particularly low and 
so the associations may not be genuine; but they 
may be genuine and so we cannot remove these 
sources from the maps. 

We can use the fact that four LBGs do lie 
within 8 arcseconds of a SCUBA source to esti- 
mate how many real SCUBA-LBG associations we 
may have missed by discarding all SCUBA sources 
at a greater distance from an LBG. On the as- 
sumption that these four are genuine associations, 
we would expect, from the results of the simula- 
tions, that 0.05-0. Ix 4 associations to have an off- 
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set of > 8 arcscconds. Since some of these four as- 
sociations are undoubtedly due to random chance, 
this gives an upper hmit of 0.2-0.4 missed identi- 
fications. 

After removing all the sources except the above, 
we measured the noise-weighted mean of the flux 
at the positions of LBGs, treating these three sam- 
ples of LBGs separately. The results are sum- 
marised in Table 2. For each sample, the mean flux 
is above zero, but clearly not at a very convincing 
level. The significance of these measurements can 
be estimated using a simple KS-test which deter- 
mines the level at which the distribution of LEG 
fluxes is inconsistent with being drawn randomly 
from the distribution of pixel values in the sub- 
mm map. The KS probabilities are given in Ta- 
ble 2. There are detections at approximately 2a- 
(KS probability -90%) for the two CFDF sam- 
ples while the Steidel LEG sample is consistent 
with zero mean flux. The distributions of sub-mm 
fluxes values for the three LEG samples are shown 
in Figure 1. 

We have also investigated the effect on our anal- 
ysis if we do not remove any sources from the 
maps. In this case, we obtain the following re- 
sults for the three samples: -0.171 ± 0.263 mJy for 
CFDF-03, 0.760 ± 0.206 mJy for CFDF-14, and 
-0.021 ± 0.110 mJy for Steidel-14. For the 14-hour 
CFDF sample, we now have a highly significant re- 
sult while for the other two fields the mean flux has 
decreased. It may seem counter-intuitive that the 
coadded flux could actually decrease if the bright 
sub-mm sources are not removed but this is sim- 
ply a result of the beam profile. The chopping 
technique creates two negative sources for every 
positive source (offset by 30 arcsecs in each direc- 
tion of RA) and so if these are not removed they 
can lead to a decrease in the coadded flux (espe- 
cially if the LBGs are clustered around the bright 
sub-mm sources). Of course, not removing bright 
sources can also lead to an increase in the coadded 
flux, as in the case for the 14-hour CFDF sample, 
if a significant number of the LBGs are close to 
the peak of the positive beam (although still with 
a > 8 arcsecs offset). 
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Fig. 1(a). — Sub-mm flux distribution for the 
Lyman-break galaxies in the 3-hour sample. Over- 
laid is the distribution of all the pixel values in 
the 3-hour submillimeter map (with all sources re- 
moved) to illustrate the fiux levels at which the 
LEG flux distribution deviates from that of the 
sub-mm map. The weighted mean of the sample 
is also shown (vertical dashed line). The sub-mm 
3-hour field had poorer weather than the 14-hour 
field and is therefore noisier. The 2 measurements 
at >5 mJy and 1 at < -5 mJy are in extremely 
noisy regions of the image (with S/N < 3a) and 
as such have a very low weighting factor in the 
mean. 



3.2. The Lyman- Br ecik GalcLxy Contribu- 
tion to the Submillimeter Background 

The contribution from LEGs to the FIR/sub- 
mm background is still an open question and one 
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Table 1 

Lyman-Break Galaxies with Ssso/iim > 3 mJy 



Source Name 


Sssomto (mJy) 


LEG Name 


P probabilities 


CFRS14-6 


4.2 


CFDF 64601 


0.019 


CFRS14-7 


3.2 


Steidel West2-MD13 


0.096 


CFRS14-8 


3.4 


Stoidol MMD75 


0.081 


CFRS14-10 


3.5 


Steidel MMD63 


0.072 



Table 2 

SUBMILLIMETER FLUX OF LYMAN-BREAK GAL.YXIES 



LEG Sample 


Mean Flux 
(sub-mm soiirccs removed) 


KS significance 


Mean Flux 
(no sources removed) 


CFRS03+00 CFDF 


0.414 ± 0.263 


94.4 % 


-0.171 ± 0.263 


CFRS14+52 CFDF 


0.382 ± 0.206 


88.7 % 


0.760 ± 0.206 


CFRS 14+52 Steidel 


0.0141 ± 0.110 


18.4 % 


-0.021 ±0.110 




Fig. 1(b). — Same as Fig. 1(a) but for the Lyman- 

brcak galaxies in the 14-hour CFDF sample. Over- 
laid is the distribution of all the pixel values in the 
14-hour submillimeter map (with all sources but 
the one source in Table 2 removed). The weighted 
mean (vertical line) is also shown. 



which can be addressed by these data. The bright 
sub-mm selected sources {Ss50ij.m > 3 mJy) pro- 
duce 20-30% of the background energy at 850/im 
(Barger et al. 1998, 1999; Hughes et al. 1998; 
Blain et al. 1999; Eales et al. 2000; Cowie Barger; 
Small et al. 2002), and therefore if the LBGs are 
the same population as the Ss5o^m < 3 mJy ob- 
jects they could be responsible for the remaining 
70-80%. Adelberger & Steidel (2000) estimate 
that [/y-selected galaxies with mieoo < 27 could 
easily have produced ~75% of the sub-mm back- 
ground at 850/im. However, in this picture the 
bulk of this energy is emitted by LBGs too faint 
to be included in our optical sample and at sub- 
mm fluxes too faint to be statistically detected 
using our method. Chapman et al. (2000), who 
targeted brighter LBGs (7^ <24.5) which are ex- 
pected to be stronger sub-mm emitters, concluded 
that these objects produce a negligible contribu- 
tion to the background. 

Given that none of our three LBG samples has 
been strongly detected as a population at 850/Ltm 
we estimate upper-limits to the background con- 
tribution that these galaxies make at this wave- 
length. Given the la standard error of the coad- 
ded flux measurements in Table 1 we can estimate 
a 3 o" upper limit for the background contribution. 
For the CFDF samples, in the redshift range 2.9 
< ^; < 3.5 the 3-hour field contributes < 3.2% to 
the background and the 14-hour field contributes 
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Fig. 1(c). — Same as Fig. 1(a) but for the 
Lyman-break galaxies in the 14-hour Steidel sam- 
ple. Overlaid is the distribution of all the pixel 
values in the 14-hour submillimctcr map (with all 
sources but the three listed in Table 2 removed). 
For consistency the flux for each object was taken 
to be the flux at the position of the LBG, including 
the three LBGs possibly associated with SCUBA 
objects at > Su. These three objects have sub- 
mm fluxes (measured at their sub-mm location) 
of < 3.5 mJy and therefore because of the offset 
between the LBG and SCUBA positions have re- 
covered fluxes of ~2 mJy. 



< 2.8%. For the Steidel et al sample, which is 
selected over a larger redshift range of 2.4 < 2 < 
3.4, the contribution is < 5.1%. 

However, because of the negative iiT-correction, 
the observed 850/im flux is constant with redshift 
(for approximately 0.5 < 2; < 6, depending on 
cosmology) and the contribution of LBGs to the 
submillimeter background is coming from a wide 
range of redshifts, not just the tight range over 
which the Lyman-break technique selects galax- 
ies. To determine the total contribution we must 
include the sub-mm flux from LBGs outside the 
selected redshift range. To do this we assume that 
LBGs have a constant comoving number density 
between redshifts 1 and 5, and then, by integrat- 
ing the comoving volume increment (dVc/dz) over 
this range, we can scale the background contribu- 
tion accordingly. Thus, upper limits on the total 
contribution to the background at 850/im becomes 

< 20% for CFDF-03, < 16% for CFDF-14 and < 
19% for the deeper Steidel et al sample in the 14- 
hour field. 

3.3. Lyman Break Galcixies Identified 
with SCUBA Sources 

The optical and near-IR counterparts to the > 
3 mJy SCUBA sources are discussed in detail in 
Webb et al. (2001b, in preparation) and Clements 
et al. (2001, in preparation), along with the iden- 
tification process. There are four SCUBA sources 
whose best identification is a Lyman-break galaxy, 
three in Steidel et al's sample and one from CFDF- 
14. 

It is interesting to note that the three SCUBA 
sources idcmtified with Steidel et al objects are 
part of a chain-like structure in the sub-mm map 
and are all within 40 arcseconds of each other. The 
remaining source, associated with an object in the 
CFDF-14 sample (but not present in the Steidel et 
al sample) is also found in the same general area 
approximately 1 arcmin to the west. 

West2-MD13, which is identified with 14.7, is 
the only object of the three possible Steidel LBG 
identifications with a spectroscopic redshift. This 
object is a quasar which lies in the largest single 
over-density within the redshift range to which the 
Steidel et al survey was sensitive. There is rea- 
son to expect that SCUBA sources might be as- 
sociated with over-dense regions in the early uni- 
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verse, (see Discussion). For example, Chapman 
at al. (2001b) observed an over-density of LBGs 
at ^=3.09 (Steidel et al. 1998) with SCUBA and 
detected a correspondingly high surface density of 
sub-mm sources. 

4. Dust Properties of Lyman-Break Galax- 
ies 

We have obtained an upper limit {2a) of ~0.4 
mJy on the average 850/xm flux of LBGs for both 
CFDF samples. We converted this to an upper 
limit on the dust mass using the following formula 
(Hildebrand 1983): 

= 7 T^T (Ij 

l^d\Vern)ij\Vem-, J- d) 

where {l-\-z)D is the luminosity-distance and Sobs 

is the flux. We assumed that all the LBGs are 
at a redshift of 3, which means that the emitted 
frequency, Vem, is 1-41 x 10^^ Hz, equivalent to 
212.5/im. We assumed that the dust mass opacity 
coefficient, Kd{v), has a value of 0.077 m^ kg"-*^ 
at 850/xm and extrapolated this to the shorter 
wavelength using the formula Kdiv) oc with 
a dust-emissivity index, (3, of 2, for which there 
is now strong evidence (Dunne & Eales (2001) 
and references therein). To compare this upper 
limit with the dust masses of nearby galaxies, we 
also calculated dust masses for the 104 galaxies 
in the SCUBA Local Universe and Galaxy Sur- 
vey (SLUGS, Dunne et al. (2000)) using the same 
formula. This survey contains a variety of galaxy 
types, from ULIRGs like Arp 220 to galaxies which 
are much more typical of the normal galaxy pop- 
ulation. To avoid the uncertainty in /? making 
the comparison between the high and low-redshift 
galaxies dubious, we used the multi-wavelength 
data that exists for the SLUGS galaxies to esti- 
mate the fluxes of the galaxies at the same fre- 
quency used for the LBGs. Therefore, even if the 
value of (3 we have assumed to calculate the dust 
mass opacity coefflcient is incorrect, the relative 
dust masses of the LB and SLUGS galaxies will 
be correct. The value of the dust mass opacity 
coefficient at any frequency is still quite uncer- 
tain (Alton et al. 2000), and so our absolute dust 
masses may not be correct, but again our compar- 
ison of the dust masses of low- and high-redshift 
galaxies will be valid, as long as the properties of 



dust are not radically different at low- and high 
redshift (Eales & Edmunds 1996). 

Since the emitted frequency for the LBGs is 
fairly close to the peak of the dust spectral energy 
distribution, our dust mass estimates are much 
more sensitive to the assumed dust temperature 
than at lower frequencies. Dunne & Eales (2001) 
have recently shown that, for the SLUGS galax- 
ies, the mass- weighted temperature of the dust 
(the correct one to use in the formula above) is 
relatively low, even for extreme ULIRG's like Arp 
220. They derive an average mass- weighted tem- 
perature for the sample of21.3±0.5 K. We have 
used this value to derive the dust masses for the 
LBGs and the individual values of mass-weighted 
temperature for the SLUGS galaxies. Fig. 2 shows 
the results for the three different cosmologies. Al- 
though the mass limit for the LBGs depends quite 
strongly on the cosmological model, it is clear that 
their dust masses must be no larger than those 
of nearby galaxies. The significance of these dust 
masses will be discussed by Dunne et al. (in prepa- 
ration). 

5. The Cross- Correlation Function Be- 
tween SCUBA Sources and Lymein- 
Break Galetxies 

We measured the SCUBA-LBG angular cross- 
correlation function in each field with a separate 
analysis for each sample in the 14-hour field. Fig- 
ure 3 shows the positions of the LBGs and SCUBA 
sources for all three LBG samples. First, we cal- 
culated the angular separations of all the A^^ x Ni 
SCUBA-LBG pairs in each field and then divided 
these into bins of angular distance. This procedure 
gives the basic function from the data, Nsi{9i) or 
LS. 

We then generated a list of random galaxy posi- 
tions for both the SCUBA sources and LBG galax- 
ies taking account of sensitivity variations in the 
images. Since we have two different sets of sources 
- SCUBA sources and LBGs - it is necessary to 
take account of the selection effects for both sets 
of sources. The Lyman-break catalogues are sub- 
sets of larger catalogues drawn from a much larger 
area of sky than the SCUBA fields, and we gen- 
erated lists of 5000 random positions for Lyman- 
break galaxies on the assumption that the sen- 
sitivity of all the catalogues is constant over the 
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SCUBA fields. 

The sensitivity of each SCUBA field is, how- 
ever, very inhomogeneous (Eales et al. 2000; Webb 
et al. 2002). To generate lists of random positions 
for the SCUBA sources that take account of the 
variation in sensitivity of each SCUBA image, we 
adopted the following procedure. Using the best- 
fit submillimeter source counts from our survey 
(Eales et al. 2000; Webb et al. 2002) we produced 
lists of artificial random SCUBA sources. For each 
source, we produced a random position on the as- 
sumption that there is no variation in sensitivity 
across the fields. Wc then used the noise images 
for each field (Eales et al. 2000; Webb et al. 2002) 
to determine whether each artificial source would 
have been detected in our survey. In this way we 
produced a list of 5000 artificial SCUBA sources. 

We calculated the angular separations of the 
real SCUBA sources and the Nr, artificial LBGs, 
giving the function SRl ■ We did likewise with the 
real Lyman-break galaxies and artificial SCUBA 
sources (LRs), and finally the artificial SCUBA 
sources and artificial Lyman-break galaxies, giving 
the function RlRs- These functions must all be 
normalized to the same number of pairs as used to 
calculate LS. 

We used two possible estimates of w{0). The 
first of these is the Landy & Szalay formalism 
(Landy & Szalay 1993): 



Fig. 2. — Histogram of the dust masses of galax- 
ies from the SCUBA Local Universe and Galaxy 
Survey (Dunne et al. 2000), superimposed with 
vertical dashed lines showing the upper limit on 
dust mass for the LB galaxies obtained in this pa- 
per. The three lines correspond to different cos- 
mologies: from left to right, (Qq = 1, = 0), 
{Qo = 0.2, A = 0.8), {Qo = 0, A = 0). 



... LS-iLRs + SRL)+RLRs 
"^') = Rjfs 



and the second is the Hamilton formula (Hamilton 
1993) 



LS X RlRs , 

W(0) = -r-::: 1 



(3) 



LRs xSRl 
Both functions gave comparable results. 

The final complication is that of the "integral 
constraint". If w{9) is estimated from an image, 
the integral 

■^J J West{0)dQidn2 (4) 

will necessarily be zero, even though the same in- 
tegral of the true correlation function will not be 
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zero for any realistic image size (Groth & Peebles 
1977). As in Roche & Eales (1999), we assumed 
that the observed angular correlation function is 
given by: 

w{e) = - c) (5) 

C can then be calculated from 



4 



-0.8 



Nr.Nr 



(6) 



in which 0ij is the angular distance between the 
i'th artificial Lyman-break galaxy and the j'th ar- 
tificial SCUBA source. For the 14-hour field this 
was calculated to be 0.0174 and for the 3-hour 
field, 0.0158. 

Fig. 4 shows our estimates of w(^) for both 
fields. To estimate our errors we used the boot- 
strap method as outlined by Barrow, Bhavsar, & 
Sonada (1984). We then fitted Equation 5 using 
these values of C with the uj{0) results from each 
field and determined a cross-clustering amplitude 
of 5.2 ± 2.9 arcsec° for the Steidel 14-hour sam- 
ple, 1.1 ± 4.4 arcsccO *^ for CFDF-14 and and 2.3 
± 3.8 arcsec°'* for CFDF-3, where the errors are 
simply errors in the chi-squared fit. At separations 
of 10 arcsecs these amplitudes correspond to (jj{6) 
= 0.82 ± 0.46 for the Steidel sample, u>{e) = 0.17 
± 0.70 for CFDF-14 and w((9) = 0.36 ± 0.60 for 
CFDF-3. Clearly, the result from the Steidel et al 
sample in the 14-hour field, which is the largest 
sample of the three, is the most (only) significant 
result, and is quite large. 

In order to investigate the possibility that the 
14-hour field is an unusually clustered region we 
calculated the auto-correlation function of the 
Steidel et al LBGs. The result is consistent with 
the results found in Giavalisco et al. (1998) using 
many fields over a larger area. For the brighter 
CFDF sample Foucaud et al. (in preparation) 
calculated an amplitude of 4.0 ± 0.7 arcsec^'^ for 
CFDF-14 which is slightly lower than the aver- 
age of the three CFDF fields (CFDF-14, CFDF-3 
and CFDF-22) of 4.7 ± 1.2 arcsec°-^. Therefore, 
the 14-hour field docs not appear to be excep- 
tional. The amplitude of the SCUBA-LBG cross- 
correlation function for CFDF-3 is comparable to 
the amplitude of the auto-correlation function of 
the brighter LBGs in the CFDF sample. 

This analysis was performed assuming that the 
LBGs and SCUBA sources were not the same ob- 
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Fig. 3(a). This figure shows the relative po- 
sitions of the SCUBA sources and LBGs in the 
CUDSS-3 field. Diamonds correspond to LBG's 
and crosses to SCUBA sources. The sensitivity 
of the SCUBA map varies substantially over the 
entire field while the LBGs are considered to be 
uniformly selected. 
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Fig. 3(b). — Same as figure 3(a) but for the 
CUDSS-14 field and the CFDF 14-hour field LEG 
sample. 
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ject and so those objects within 8 arcsccs of cacli 
other were included in the correlation signal. To 
be certain that these few objects were not biasing 
the clustering signal we redid the analysis, with 
the added restriction of0 > 8 arcsecs. We find the 
following amplitudes: for the Steidel sample, 6.7 
± 1.9 arcscc"•^ for CFDF-14, 0.4 ± 3.4 arcsec"-^ 
and for CFDF-03, 3.8 ± 3.9 arcsec°-*. It does not 
appear, therefore, that the possible LBG-SCUBA 
associations are strongly affecting the correlation 
results. 

In this analysis we have used the complete 
CUDSS SCUBA catalogue of 50 sources. As 
this sample includes all sources detected above 3a 
one might worry that a large number of spurious 
sources are contaminating the clustering analysis. 
However, in Eales et al. (2000) and Webb et al. 
(2002) we discuss the number of expected spuri- 
ous sources in the sample, using Gaussian statis- 
tical arguments and arguments based on an anal- 
ysis of the noise in the SCUBA maps. We have 
concluded that about 2-3 sources in each field are 
spurious, or a total of 4-6 sources in the combined 
sample (or approximately 10%). At this level these 
sources are not expected to significantly alter the 
cross-clustering analysis. 

However, to check this, we have repeated the 
clustering measurement using a subset of the 26 
SCUBA sources in our catalogue which were de- 
tected at > 3.5(T. For a larger sample size one 
might expect the clustering amplitude to increase 
with the removal of the spurious sources in the 
sample, however, in doing so we are substantially 
decreasing our sample size (as we must remove 
all < 3.5(7 sources, not just the 10% which are 
spurious). We measure the following clustering 
amplitudes: A=2.7±5.2 arcsec"'^ for the 3-hour 
field, A=1.2±5.4 arcsec" *^ for the 14-hour field 
and CFDF-14, and ^=5.1±5.5 arcsec°-^ for the 
14-hour field and the Steidel et al. sample, where 
the errors have again been estimated using the 
bootstrap method. The measured values are es- 
sentially the same as the amplitudes measured for 
the entire >3.0(T sample, though the uncertain- 
ties have increased substantially because of the 
decrease in sample size. Hence, it does not ap- 
pear that the analysis has been contaminated by 
spurious sources. 



6. Discussion 

6.1. Submillimeter Flux of Lyman-Break 
Galaxies 

We detect, at a level of ^2(t, a mean sub-mm 
fiux from the two CFDF samples of LBGs but we 
do not detect any flux from the Steidel sample. 
RecaUing that the CFDF 14-hour LBGs are an 
optically bright subset of the Steidel LBGs this 
could be taken as evidence of a positive correlation 
between observed optical and sub-mm flux. How- 
ever, a plot of /-band magnitude versus recovered 
sub-mm fiux shows no such correlation (Figure 5). 
Similarly, no correlation is present when TZ-hand 
magnitude for the Steidel sample is plotted versus 
sub-mm flux. In addition, there is no statistical 
detection of the sub-mm fiux of a brighter sub- 
sample (7?. < 25 ) of the Steidel et al list (which 
agrees with the CFDF-14 catalogue at approxi- 
mately the 70% level). 

On closer inspection we find that the CFDF- 
14 measurement is completely dominated by the 
detection of one LBG and when this object is re- 
moved from the list a positive mean flux is de- 
tected only at the la level. Although there are 
three detections of LBG's at >3a in the Steidel 
list they do not lift the mean flux of the entire 86 
galaxy sample to a significant level. In CUDSS-3, 
where no LBG was detected above 3a the positive 
detection of fiux is not due to any single object. 

It is not entirely surprising that we see no corre- 
lation between the sub-mm and observed optical 
fiux. Many authors (Shapley et al. 2001; Adel- 
berger & Steidel 2000) have claimed that galaxies 
with higher intrinsic UV luminosities (after cor- 
recting for reddening) contain larger dust masses. 
However, this does not necessarily translate into a 
correlation between observed optical and sub-mm 
fiux. Indeed, as with our results, Shapley et al. 
(2001) found no correlation between dust content 
and observed optical fiux. 

Two other groups have attempted to directly 
measure the sub-mm fiux of high-redshift star 
forming galaxies. Chapman et al. (2000) carried 
out a targeted study of 8 Lyman-break galaxies, 
selected to have high UV-devived star formation 
rates and they obtain a result similar to our own. 
Although one object was detected at >3a, no fiux 
was detected statistically from the sample as a 
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whole once this object was excluded. However, 
in a recent conference proceeding (Chapman et al. 
2001a), with a larger sample of 33 red and high- 
star forming LBGs, they report a statistical de- 
tection of S85o=0.6 ±0.2 mJy with marginal de- 
tections for several individual red LBGs. 

Peacock et al. (2000) followed a similar ap- 
proach to this paper for starburst galaxies in the 
Hubble Deep Field and detected a mean flux of 
0.2 ± 0.04 mJy (for a SFR of 1 /i-^Moyr"!) with 
the mean flux increasing with SFR. To directly 
compare this result with our own we must first 
convert our statistical measurement to units of 
submillimeter flux per unit star formation rate. 
To determine the UV-estimated star formation 
rates for the LBGs in our three samples we fol- 
low the method outlined in Peacock et al. (2000) 
and Madau et al. (1996). Approximately 60% of 
the Steidel et al. LBGs in our field have spectro- 
scopic redshifts. For those without spectroscopic 
redshifts, and for the CFDF LBGs, we assumed a 
redshift of z=3. We chose a flat, f2A=0.7 cosmol- 
ogy. We flnd the following: 



SFR 



for the Steidel et al. 14-hr sample 

5850^m/mJy = 0.015 ± 0-022 
for the CFDF 14-hr sample 



T (9) 



<S'850/.m/mJy = 0.065 ± 0.034 



for the CFDF 3-hr sample 



5850^m/mJy = 0.13 ±0.06 



SFR 



/i-2M0yr- 



T (10) 



SFR 



/i-^Mgyr- 



T (11) 



Of our three LBG samples only the statistical 
measurement from CFDF-03 is consistent with the 
results of Peacock et al. (2000). The Steidel et al. 
sample and the CFDF-14 sample both have much 
lower detections. Peacock et al. (2000) converted 
the results of Chapman et al. (2000) to these units 
and found: 



^850Mm/mJy = 0.13 ±0.14- 



SFR 



— 

'/i-2MQyr-i ^ ' 
which is consistent with the Peacock et al. (2000) 



result, although also consistent with no detection 
of submillimeter flux. It appears that most LBGs 
are not strong sub- mm emitters, that is with 0.3 
< Ss5Q^rn < 3 mJy, and therefore the strength of 
a statistical submillimeter detection may depend 
strongly on the properties of the specific LBG sam- 
ple observed. 

In figure 5 we plot the mean 850/im fiux as a 
function of Z/F-estimated star formation rate for 
all three LBG samples. Though, to the eye, some- 
what suggestive of a rise, these data are consis- 
tent with no correlation between the star forma- 
tion rate and the submillimeter flux. 

As discussed in §4 there are four Lyman- 
break galaxies which are possibly identified with 
a SCUBA source at >3cr, one from the CFDF 
14-hour sample and three from the Steidel ct al 
sample. There are very few Lyman-break or high- 
redshift star forming galaxies (from all studies) 
which are known to be sub-mm bright and there- 
fore conclusions regarding any unifying properties 
are difficult. There is some suggestion however 
that these objects have extremely red colors com- 
pared to the average for the population (Chapman 
et al. 2000, 2001a). As discussed in Shapley et al. 
(2001) the more intrinsically luminous LBGs ap- 
pear to be dustier, with redder optical colours, 
and should therefore be brighter in the sub-mm. 
However, we see no such trend with the four LBGs 
identified with bright SCUBA objects in our sam- 
ple although we have very limited photometric 
information, particularly in the near-infrared, on 
these faint objects. 

We estimate that the LBGs from the two CFDF 
samples could be producing at most 20% of the 
background at 850/xm, when integrated over all 
redshifts. The Steidel LBGs, which are a deeper 
sample have an upper limit to their contribution 
of 19%. However, as outlined in Chapman et al. 
(2000), the LBGs that are most likely do domi- 
nate the sub-mm background, the highly reddened 
galaxies, are less likely to be detected in these op- 
tical LBG surveys and therefore our sample may 
be biased to sub-mm faint objects. 

6.2. The SCUBA-LBG Cross-Correlation 
Function 

The angular cross-correlation functions be- 
tween the SCUBA sources and the three LBG 
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samples arc presented in Figure 4. There arc some 
interesting results from this clustering analysis. 
The first is that, though we measure consistent 
clustering amplitudes for all three samples, our 
strongest LBG-SCUBA cross-clustering signal is 
detected in the Steidel et al. sample in the 14-hour 
field. Recent results from Giavalisco & Dickinson 
(2001) and Foucaud et al. (in preparation) have 
shown clustering segregation with UV luminosity 
for the auto-correlation function of the LBGs, such 
that the brighter objects are more strongly clus- 
tered. Thus one might posit that we should see 
a stronger clustering signal in the CFDF samples. 
Unfortunately, the smaller numbers in the brighter 
CFDF samples mean that even if a stronger signal 
was present it would be harder to detect. 

A second result is that the measured ampli- 
tude of the cross-clustering between the 14-hour 
SCUBA sources and the Steidel LBGs is larger 
than for the auto-clustering of the LBGs them- 
selves, though certainly consistent within the er- 
ror. The strength of the LBG-SCUBA angular 
cross-clustering is essentially a lower limit on the 
true spatial clustering as it is projected over a 
broad redshift range. The actual spatial cross- 
clustering is expected to be even higher, since the 
range in the redshift distribution of the SCUBA 
sources is much broader than that of the LBGs. As 
a further check we repeated the cross-correlation 
analysis after removing SCUBA sources with se- 
cure IDs with z < 2 (Eales et al. (2000), Webb 
et al, in preparation, Clements et al, in prepara- 
tion) and found the correlation signal increased as 
would be expected for real clustering. 

By assuming a redshift distribution for the 
SCUBA sources (including those with z < 2) we 
can estimate ro for the spatial LBG-SCUBA cross- 
correlation function following the procedure in (for 
example) Efstathiou et al. (1991). Although the 
redshift distribution is still highly uncertain, the 
results of many groups are not inconsistent with 
a median redshift near z=3 and with only a small 
fraction of objects below z=2 (see the review by 
Dunlop (2001) and references therein). We there- 
fore take our general redshift distribution to be a 
Gaussian centered at z=3 with (t=G.8. We adopt 
the redshift distribution of Giavalisco et al. (1998) 
for the Steidel LBGs and f2m=0.3 and f2A=0.7 cos- 
mology. Given these parameters we find ro=11.5 
± 3.0 ± 3.0 h-^ Mpc for the Steidel et al sam- 



ple. The first error is simply the statistical error 
calculated from ui{0) and 6ui{0). The second er- 
ror is systematic and is estimated by varying the 
redshift distribution of the SCUBA sources from 
z=2.5 to z=3.5, and Az from 0.6 to 1.1. The 
redshift distribution of the CFDF LBG samples 
is slightly different than that of Steidel et al. We 
estimate the CFDF distribution as a Gaussian of 
the same standard deviation (0.24) but centered 
at z =3.2. For the CFDF-14 amplitude of 1.1 ± 
4.4 arcsec°-® we find ro=4.5 ± 7.0 ± 5.0 h'^ Mpc 
and for CFDF-3 (with an amplitude of 2.3 ±3.8 
arcsec°-S) we find ro=7.5 ± 7.0 ± 5.0 Ir'^ Mpc. 

The two values of Tq determined from the 
SCUBA-CFDF results are comparable (within 
their large uncertainties) to those measured for the 
LBG auto-correlation function in the two fields. 
In the CFDF, for fi^^O.S and f2A=0.7 cosmology, 
Foucaud et al. (2002, in preparation) estimated 
To = 6.4 ± 0.3 Mpc for the 3-hour field and 
ro = 5.1 ± 0.5 h-^ Mpc for the 14-hour field. 

We consider the Steidel LBG-SCUBA cross- 
correlation result to be the most secure. For Stei- 
del LBGs with TZ < 25.5, Giavalisco & Dickin- 
son (2001) found ro = 3.2 ± 0.7 for the auto- 
correlation function, for the same cosmology which 
is smaller than our result for this field, though 
within the uncertainty range. One simple argu- 
ment suggests that a higher value for the ampli- 
tude of the LBG-SCUBA cross-correlation func- 
tion than for the Lyman-break galaxies themselves 
would not be unexpected. The high value of the 
auto-correlation function for Lyman-break galax- 
ies alone has been explained by the large values of 
bias expected for rare systems in the early universe 
(Kaiser 1984; Giavalisco et al. 1998). Indeed, the 
more luminous LBGs have been shown to be more 
highly clustered (Giavalisco & Dickinson 2001). 
Luminous SCUBA sources are much rarer objects 
than Lyman-break galaxies, and so the bias values 
may be even higher. This argument breaks down 
if SCUBA sources are the result of rare or short- 
lived stochastic processes in the universe, in which 
case they need not be highly-clustered at all, with 
themselves or with the LBGs. 

Webb et al. (2002) and Scott et al. (2002) 
have measured the auto-clustering strength of 
SCUBA sources. Though hampered by small ar- 
eas and numbers the results are consistent with 
strong auto-clustering, at least as strong as the 
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auto-clustcring of LBG galaxies and EROs. Webb 
at al. (2002) find an angular correlation ampli- 
tude of 4.4±2.9 6~^-^. Assuming the same redshift 
distribution as above this may bo inverted to a 
spatial amphtude of 12.8±7.0 h~^Mpc. 

It is tempting to draw an analogy with the 
universe at low redshift, where the amplitude 
of the correlation function for rare clusters of 
galaxies is ~18 times higher than that for galax- 
ies themselves, with the amplitude of the cross- 
correlation function between clusters and galax- 
ies being midway between these values (Bahcall 
1988). Circumstantial evidence in favour of this 
idea is the discovery of clusters of submillimetre 
sources around the extremely rare high-redshift ra- 
dio galaxies (Ivison et al. 2000), suggesting again 
that SCUBA sources preferentially form in much 
rarer environments than the more Lyman-break 
galaxies. However, from the results of this work 
we may only say that the cross-clustering between 
the SCUBA sources and LBGs is at least consis- 
tent with the strengths of the self-clustering in the 
individual populations. 

A cross-clustering signal has also been detected 
between very bright SCUBA objects and Chandra 
sources (Almaini et al. 2001) with an even larger 
amplitude than found in this work. X-ray bright 
objects, as with sub-mm bright objects are rela- 
tively rare and therefore, following the same logic 
as above should also be very highly clustered. 

The possible angular cross-clustering signal 
measured in the Stcidcl 14-hour field is sugges- 
tive evidence that many of the SCUBA sources in 
this field are indeed at the same high-redshifts as 
the LBGs. The smaller angular cross-clustering 
signal in the 3-hour field may indicate that the 
SCUBA sources in this field do not lie at the same 
rcdshifts as the LBGs and would suggest a differ- 
ence in the redshift distribution of SCUBA sources 
in the 3 and 14-hour fields. 

The positive clustering signal should be remem- 
bered when attempting to determine identifica- 
tions for SCUBA sources by positional coincidence 
or when observing individual LBGs with a large- 
beam telescope such as the JCMT. Positive sub- 
mm flux could erroneously be associated with a 
near-by LBG rather than the object actually pro- 
ducing the emission which might be undetected 
at optical wavelengths. Indeed, this may be the 
case with our own LBG identifications and those of 



other authors, in particular sub-mm bright LBGs 
which do not show unusual colours or luminosi- 
ties. 

7. Conclusions 

We have used the 850 /xm maps from the 
Canada-UK Deep Sub-mm Survey to study (i) the 
sub-mm flux and dust properties of Lyman-break 
galaxies and (ii) the angular correlation between 
Lyman-break galaxies and SCUBA sources. We 
obtain the following results: 

1. We marginally detect (at the 2a level) sub- 
mm flux from Lyman-break galaxies in the CFDF- 
14 and CFDF-03 samples but wc do not detect flux 
from the Steidel et al sample. The flux levels are: 
0.382 ± 0.206 mJy for the 14-hour field and 0.414 
± 0.263 mJy for the 3-hour field. Further, wc show 
that possibly because of LBG-SCUBA clustering 
SCUBA sources not identified with a LBG galaxy 
must be removed from the map before a proper 
analysis can be performed. 

2. Lyman-break galaxies are the best optical 
identification for four SCUBA sources although it 
is possible that some of these identifications may 
be incorrect. There are indications that these ob- 
jects may lie in a region of spatial over-density. 

3. An upper limit for the dust mass of Lyman- 
break galaxies was calculated from their sub-mm 
fiux results and we conclude that these masses can 
be no larger that those of near-by galaxies. 

4. The SCUBA-LBG correlation function was 
measured for all three sample of Lyman-break 
galaxies. Wc found a high-amplitude ro= 11.5 ± 
3.0 ± 3.0 Mpc for the Steidel et al sample, 

= 4.5 ± 7.0 ± 5.0 h-^ Mpc for CFDF-14 and To = 
7.5 ± 7.0 ± 5.0 h-^ Mpc for CFDF-03, (where the 
first error is statistical and the second systematic) . 
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Fig. 3(c). — Same as for figure 3(a) and (b) but for 
the CUDSS-14 field and the Steidel et al 14-hour 
field LBG sample. 
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Fig. 4. — Angular cross-correlation function 
(SCUBA-LBG) for the three Lyman-break sam- 
ples. The triangles are for the 3-hour field, the 
filled circles are for the Steidel et al. sample of the 
14-hour field and the open circles are for the CFDF 
sample of the 14-hour field. The upper, bold solid 
line is the LBG-SCUBA cross-correlation function 
measured in this paper (not including the inte- 
gral constraint) for the Steidel et al sample. This 
fit used only the inner four points since beyond 
about 150 arcsoconds edge effects become impor- 
tant. The middle dotted- line corresponds to the 
auto-correlation function of the bright (Sgso/jm > 
3 mJy) SCUBA sources from Webb et al. (2002). 
The lower solid-line is the auto-correlation func- 
tion for the Lyman-break galaxies from Giavalisco 
& Dickinson (2001). 
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Fig. 6(a).— The I-band flux of the CFDF LBGs 
for both the 3-hour and 14-hour fields versus their 
recovered 850 /im flux. 
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Fig. 5. — The mean 850/xm flux as a function of 
J/V'-estimatcd star formation rate. The squares 
and diamonds correspond to the CFDF-14 and 
CFDF-03 samples respectively. The triangles rep- 
resent the Stcidcl-14 list. Because of our small 
sample size we are restricted to only two bins for 
the CFDF fields and four for the Steidel field. The 
points are consistent with no detected rise in flux 
with star formation rate. 
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Fig. 6(b).— The R-band flux of the Steidel et al 
14-hour field LBGs versus their recovered 850 yum 
flux. 
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